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ABSTRACT 


A Monte Carlo computer simulation was developed to model 
hypothesized electro-optic laser communication systems 
operating in the middle ultraviolet region of the spectrum 
called the solar blind. By assuming various source, propa- 
gation, and detector characteristics as well as certain 
performance parameters it is possible to predict the effec- 


tive ranges and operating characteristics of such a system. 
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I. INTRODUCTION E 


A. HISTORY 

Historically optical communications systems were devel- 
oped over line-of-sight paths due to the propagation charac- 
teristics of light. In the United States Navy line-of-sight 
communications became somewhat neglected with the advent of 
uhf radio technology. For years the only optical communica- 
tlon systems in the navy were the 12 inch (and for a time 
24 inch) hand modulated signal lamps, yardarm blinkers, flag 
hoist, semaphore, and NANCY (a hand modulated infrared sys- 
tem). Some experimentation was conducted in the late 1940's 
and 1950's using cesium lamps operating in the near infrared 
(ir) and non-line-of-sight or over-the-horizon (OTH) commu- 
nication was demonstrated. Interest waned, however, and 
progress stood still. 

Technology has not stood still and today virtually any 
radio signal can be monitored and the transmitter localized 
using a direction finding receiver. While communication 
systems must be reliable and secure, sometimes it is impor- 
tant that they be covert (only detectable by friendly 
forces). Thus armed with a complete set of radio communica- 
tion equipment, the unit commander nevertheless may find 
his forces in EMCON, talking to each other via flashing 


light: 
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B. VISIBLE LIGHT 

In the 1960's several attempts were again made to demon- 
strate the feasibility of OTH and over-terrestrial-obstacles 
communication using atmospheric scattering [1-3]. These 
efforts utilized visible search lights, flash lamps, and 
lasers, and communication at night was successfully demon- 
strated over distances up to fifty kilometers. 

Visible radiation is highly attractive because of the 
low extinction coefficients (depending on weather condi- 
tions). Unfortunately it is seriously degraded during day- 
light hours due to the background radiation from the sun 


which tends to interfere with detection. 


C. INFRARED LIGHT 

Recently, efforts have been made to study extended line- 
of-sight communications utilizing lasers operating in the 
near ir. One such effort utilized a 1.06 y laser and 
employed remotely piloted vehicles (RPVs). Communications 
at this wavelength appeared quite promising and ranges from 
25 to 150 nautical miles have been predicted [4]. 

In this region of the spectrum, however, there is 
noticeable background noise due to scattering of incident 
solar radiation and to emission by atmospheric particles 
heated by incident radiation. It should be noted that the 
daylight background visible radiance is due primarily to 
scattered radiation and night radiance is caused mainly by 


scattered radiation and moon light. In the near ir daylight 


IRE 


background radiance is due mainly to atmospheric emission, 
while at night aurora and afterglow are the predominate 


sources. 


D. ULTRAVIOLET LIGHT 

There is a portion of the spectrum, however, that has 
been studied little and may offer much. This is the middle 
ultraviolet (uv). Atmospheric absorption in the .2 to .3 H 
region is caused primarily by ozone [5-6]. The relatively 
dense layer of ozone (at approximately 22 km) prevents vir- 
tually all of the radiation in the 220 to 280 nm region from 
reaching the surface of the earth. Thus the limiting noise 
in an electro-optical communication system operating in the 
middle uv would be self-noise rather than atmospheric noise. 
Another critical feature to be considered is the attenuation 
properties of the intervening atmosphere. In the near uv 
(wavelengths between .3 and .4 y) the primary attenuation 
mechanism at low altitude is scattering. This can be divided 
into Rayleigh scattering and particulate scattering by 
atmospheric aerosols. Rayleigh scattering, inversely pro- 
portional to the fourth power of wavelength, increases quite 
rapidly as the wavelengths become shorter. In practice, it 
has been difficult to determine the particulate scattering 
in the near uv. For clear sky conditions, however, a rea- 
sonable estimate of attenuation can be made by assuming that 
the entire attenuation is due to Rayleigh scattering [2]. 

In the middle uv particulate scattering is very sig- 


nificant. At high altitudes the large ozone concentration 
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gives rise to very large absorption coefficients. The con- 
centration of ozone generally decreases as the altitude is 
decreased. Yet even at sea level the concentration is still 
high enough to have a significant effect. These propaga- 
tion characteristics indicate that only short range commun- 
ication systems operating in the middle uv are likely to be 
useful. Because of the inherent scattering characteristics 
of middle uv radiation, along with the minimal background 
problem in the solar blind region, short range non-line-of- 
sight communications utilizing broad band light sources 
have been developed [7-8]. 

Recently a number of reasonably efficient uv lasers 
have been developed with energies of .1 to one joule per 
pulse being reported [9-10]. These represent a new genera- 
tion of radiation sources to be considered for use in the 
solar blind. A number of detectors and filters have been 
available for some time. It was these considerations that 


led to this study. 


E. PURPOSE 
The purpose of this effort is to develop an analytical 
tool for estimating the ranges for which optical communica- 


tion is possible in the middle uv. 
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II. BACKGROUND 


A. INTRODUCTION 

A typical communication system consists of a transmit- 
ter, a receiver, and the intervening space between them. 
The characteristics of the transmitter and the receiver 
are highly dependent on each other and on the properties 
of the propagating medium. Each of the components may, 
however, be treated as a distinct and relatively separate 
entity. The radiation source, or transmitter, must be 
suitably modulated to carry information and must have 
sufficient power to overcome the attenuating mechanisms in 
the medium in which it operates. The characteristics of 
the propagating medium are not well defined and in general 
must be approximated. The detector, or receiver, must have 
sufficient sensitivity to detect and demodulate the atten- 
uated signal that is embedded in background noise. Each of 
these components is present in an electro-optical communica- 
tion system and is discussed below. It will be seen that 
1) advances in sources dictate a new look at the problem; 
2) uncertainties in propagation characteristics of the 
middle uv requires careful experimentation to gather suitable 
information such that a proposed system may be modeled; and 
that 3) certain aspects of detectors, such as filters, 
Critically determine the performance parameters of such a 


system. 
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B. SOURCES 

Ultraviolet sources suitable for communication pur- 
poses include omnidirectional broad band lamps such as 
xenon or mercury flash lamps, and lasers. Experiments 
utilizing omnidirectional sources have been carried out 
for several years at the Aeronautical Research Associates 
of Princeton (ARAP), Princeton, New Jersey [7-8]. Uti- 
lizing a pulsed source at 4810 pps and an average power level 
of ten to thirty watts, voice communication was demonstrated 
over a distance of several hundred meters. With suitable 
modifications, it is estimated that ranges of up to three 
kilometers may be achieved with the current system. With 
such a source, the conversion efficiency from input power 
to output power is about 0.1 percent. 

New uv laser sources reported in the last year or so 
appear to have demonstrated orders of magnitude improvement 
in efficiency. Utilizing rare gas-halide eximer species, 
uv output with working efficiencies of up to ten percent 
are conceivable. Practical working models are already 
commercially available. Wavelengths in the range of 200 to 
400 nm have been reported, indicating that various wave- 
lengths in the middle uv may be attained with these sources. 
The power output of these devices typically runs from 0.1 
joules upwards, depending on the size of the device. The 
high degree of collimation of the output suggests that the 
geometric patterns for communication will be highly direc- 
tional. This is one aspect of the problem which is quite 


different from that observed for omnidirectional sources. 
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C. PROPAGATION 

The absorption coefficient for middle uv radiation 
varies by orders of magnitude as a function of wavelength. 
Scattering by particulates at these wavelengths is more 
severe than in the visible region of the spectrum. The 
combined absorption and scattering cross-sections give rise 
to extinction coefficients which have been measured to be 
in the range of 0.1 to five per kilometer [11]. Further, 
variations of an order of magnitude in extinction coefficient 
over a spectral range of less than fifty nanometers indicates 
extreme sensitivity to wavelength, due largely to the 
absorption bands of species such as ozone which appear in 
this region. 

At wavelengths less than 280 nm, the absorption coef- 
ficient of atmospheric ozone is sufficiently intense as to 
eliminate the large majority of solar radiation from the 
earth's surface [12]. This is largely due to an ozone 
concentration centered at 20 to 25 km above the surface of 
the earth. At sea level, the ozone concentration is typi- 
cally an order of magnitude less than at the higher alti- 
tudes. This gives rise to reduced solar background levels 
without the high extinction coefficients found at higher 
altitudes. Such behavior suggests that communication in 
the middle uv may be feasible for relatively short distances 
at sea level, with a minimum of interference from solar 


background radiation. 
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Further major considerations for middle uv propagation 
are the effects due to multiple scattering. Although the 
primary light beam may be attenuated by an extinction 
coefficient which is readily measured, the scattered radia- 
tion, especially after multiple scattering, may contribute 
considerable intensity to the radiation density at an 
observer.  Shettle and Green [12] present calculations for 
solar flux as a function of the sun's angle with respect to 
the zenith. Particularly at low angles of incidence, the 
diffuse radiation is found to be twenty to thirty orders of 
magnitude more intense than that due to the direct solar 
radiation. Such figures clearly emphasize the importance 
of properly treating the effects of multiple scattering 
when taking into account the effects of particulates in the 
atmosphere. 

In general multiple scattering calculations are cumber- 
some to handle and yield limited information. Typically 
these calculations are done on a computer using either an 
analytical [13] or a Monte Carlo [14-15] appraoch. The 
Monte Carlo calculations are relatively easier to model, 
but utilize a significant amount of computer time and only 
approximate information is obtained. Some generally useful 
results have been presented by Bucher [14-15] again indicating 
orders of magnitude greater radiation densities from the 
multiple scattering effects than from transmission of direct 
radiation. The available results from such calculations do 


not readily answer questions such as 1) what is the 
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scattering intensity for non-line-of-sight communication 
given an obstacle with an assumed angle between transmitter- 
obstacle-receiver; and 2) what is the actual radiation 
density as a function of field of view and/or position of 

an observer located off from the axis of propagation at a 


given distance from a collimated source such as a laser. 


D. DETECTION/DEMODULATION 

The information available from a pulsed signal depends 
on the pulse width, repetition rate, and signal-to-noise 
ratio. Considerations such as these are included in papers 
by Kennedy [16-17] in which communication through optical 
scattering channels is discussed. Considerations leading 
to a suitable вый CE ratio analysis are also pre- 
sented by Yarif [18] for electro-optical systems. Whether 
the information is pulsed or continuous wave greatly affects 
the reception characteristics of the signal. The background 
light levels coming through the filters utilized before the 
detector critically determine the signal-to-noise ratio. 
The characteristics of the filter and detector thus are 


crucial in determining the information available in a 


communication system. 
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III. STATEMENT OF THE PROBLEM 


One of the primary constraints on the development and 
deployment of any new system in the United States Navy is 
the limited fiscal resources available for such efforts. 
With the advent of new technology it is nevertheless impos- 
sible to utilize the technology to build a new system 
without first having a mission or demonstrated need for that 
system and second having a well-defined degree of confidence 
in the ability of that system to meet that need. 

An electro-optical communication system utilizing a 
laser operating in the middle uv is such a system. (It 
requires little or no imagination to visualize the advan- 
tages of rapid short range, possibly non-line-of-sight, 
optical communications over the methods presently employed 
by the Fleet.) While it is beyond the scope of this 
project to develop such a system, it is within the scope 
of this effort to formulate a model of such a system capable 
of predicting the various parameters of that system. 

The need for modeling is further enhanced by the limited 
amount of information available concerning the middle uv. 

By utilizing available information and by characterizing 
Source, propagation, filter and detector characteristics, it 
is possible to formulate a model that will yield order of 
magnitude results for any proposed system. In the process 


of modeling and analysis the need for critical experiments 
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which would allow the modeling efforts to proceed with 


a greater degree of accuracy is sure to be revealed. 
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IV. LASER SOURCES 


The past two or three years have seen a rapid develop- 
ment of lasers operating in the ultraviolet. The most 
notable recent advances have occurred with the rare gas- 
halide lasers in which an excited state eximer is formed. 
With species such as XeF*, the ground state is not an asso- 
ciated molecule, thus depopulation of the lower level is 
not a barrier to population inversion. Further, such rare 
gas-halide lasers have an inherently high quantum efficiency. 
Efficiencies relative to the energy deposited by an excitat- 
tion electron beam are as high as ten percent, while overall 
working efficiencies of up to one percent have been reported 
[9,10,19]. With higher overall efficiencies anticipated, 
such lasers represent an attractive source of radiation for 
middle uv applications. 

The observed wavelengths of several rare gas-halide 
lasers are presented in Table I. It should be noted that 
the laser output from these sources may be Raman shifted 
with reasonable high conversion efficiencies [20]. By 
passing the output from a KrF laser through ten atmospheres 
of hydrogen, conversion to the first stokes wavelength of 
25 percent and to the second stokes wavelength of ten percent 
was observed [20]. At eighty atmospheres pressure, more 
than fifty percent conversion efficiency to the first 


stokes wavelength has been observed. 


DAI 


TABLE I 


Characteristic Wavelengths of Several Rare 
Gas-Halide Lasers Along with Some Stokes 
Shifted Wavelengths Observed with Hydrogen Gas 


First Stokes Second Stokes 
Eximer Wavelength (nm) Wavelength (nm) Wavelength (nm) 
KrF 249 279 318 
XeF 351 
ΧΕΟΙ 308 
XeBr 282 322 


The practical performance of rare gas-halide lasers is 
illustrated by the commercial availability of a model which 
will produce 0.1 joule pulses at 15 pps. One table top 
model is reported to produce 1.5 joule pulses with a 125 
nanosecond pulse width. The development of these lasers is 
in its infancy. Sources with higher repetition rates and 
greater reliability are to be expected. 

Due to the rapidly expanding availability of middle uv 
lasers sources, no single presently available source is 
assumed, but rather the effort in the succeeding sections 
considers the possibility of utilizing a middle uv source 
at representative wavelengths and with an output power 
consisting of millijoule to joule pulses. Thus the behavior 


of any unforseen laser or coherent beam may be predicted. 
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V. PROPAGATION PROPERTIES OF THE ATMOSPHERE 
AND THE MULTIPLE SCATTERING PROBLEM 
A. INTRODUCTION 

A laser signal propagating through the atmosphere is 
subjected to various attenuation mechanisms. The time 
spread of a pulse can be broadened due to scattering, while 
the amplitude can be decreased due to absorption and scat- 
tering by atmospheric constituents. The inherent direction- 
ality of a laser source indicates that the position of the 
detector relative to the source should also affect the 
reception of the signal. Since a large fraction of the total 
extinction process is due to scattering, multiple scattering 
should be an important consideration affecting the propaga- 
tion characteristics for distances greater than one or two 
extinction lengths. 

A complete description of all of the effects of the 
atmosphere including inhomogeneities (clouds, terrestrial 
objects, surface effects, etc.) is beyond the scope of this 
paper. A relatively concise treatment of the atmosphere 
and the multiple scattering problem entails characterizing 
the atmosphere, simulating atmospheric effects, and tabulating 


the results of the simulation. 


B. CHARACTERISTICS OF THE ATMOSPHERE 
1. Attenuation 
In order to characterize the propagation of optical 


information in the atmosphere it was first necesary to 
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determine the extinction, absorption, and scattering coef- 
ficients. Table II was constructed based primarily on the 
information reported by Dunkelman [21]. The Rayleigh 
Scattering coefficients were adapted from Penndorf [22]. 
For the purpose of this investigation the following 


definitions and relationships apply: 


À is the wavelength 

Κλ 15 the absorption coefficient 

Kon is the Rayleigh scattering coefficient 

Κορ is the particulate scattering coefficient 

Kc is the total scattering coefficient 

K is the total extinction coefficient 

Po IE EE 

K = Kç + K, 

Ro is the ratio of Κς to K or (K/K) 

Rp is the ratio of Kop to Ko or (Kop/ Ko) 
From Table II it can be seen that Kar Κος’ and K 


tend to decrease as wavelength increases in the middle uv. 


R and R 


and consequently K S” p 


The relationship of Κορ ο” 


to wavelength is not so easily described. 
2.  Angular Distribution of ScHtteredgbEhOeNS 
Complete characterization of the propagation of 
optical information required the angular distribution or 
phase function, P(9), of the scattered photons for the 
wavelengths of interest. Figure 1 [23] was the primary 


source for phase function information. 
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FIGURE 1. Absolute Phase Function vs. Scattering Angle 
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In order to mathematically model this information, 


the Henyey-Greenstein function, Eq. l, and some of its 


modifications were studied. The Neer-Sandri function [7], 


Eq. 2, was a logical choice; however, the function used by 


Zachor [24], 


Eq. 3, was chosen for simplicity in modifying 


the author's initial phase function to properly account for 


backscattering. 


P(9,9) 


r (omg) 
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Fltting the data reported by Neer [23] with Eq. 3 


resulted in Figures 2 and 3. These figures represent the 


absolute and normalized phase functions, while Table III 


lists the coefficients, utilized for this simulation. 
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FIGURE 3. Normalized Phase Function vs. Scattering Angle 
Utilizing a Modified Henyey-Greenstein Function 
to Approximate the Data by Neer, et al. [23] 
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C. DESCRIPTION OF THE CHARACTERISTICS OF THE SIMULATION 

The purpose of this simulation was to provide approxi- 
mate results for the multiple scattering problem. By 
assuming various source, propagation, and detector charac- 
teristics, a user could in principle predict the effective- 
ness of a proposed system. 

The Monte Carlo method described in Appendix A was employed 
to generate information detailing characteristic patterns of 
photon scattering as a function of wavelength, position of 
an observer, pulse spreading, and field of view cf the 
detector being employed. The observer may be located at any 
arbitrary angle off from the axis of propagation and at any 
arbitrary distance from the source (measured in extinction 
lengths). 

Of primary interest were the characteristic patterns of 
multiply scattered radiation. The simulation provided for 
obtaining diagrams of photon flux at a given position as 
viewed by either the source or the receiver. In addition, 
the patterns received by the detector were capable of being 
limited by field of vidt restrictions. By measuring the 
photon flux at a given position, contours of equal photon 
flux were generated. These contours provided a simple 
means of visualizing the radiation patterns. 

Detector field of view considerations indicate the 
field of view required to receive a given threshold of 
photons at a fixed position as well as the maximum range at 


a given angle that a threshold number of photons can be 
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detected for a fixed field of view. Finally, pulse spreading 
information due to the time delay caused by multiple 


scattering was also available. 


D. CHARACTERISTIC PATTERNS OF MULTIPLY SCATTERED RADIATION 

Figures 4 and 5 represent the pattern of photons 
received by a 180 degree field of view detector located 
zero to nine and 36 to 45 degrees off from the axis of 
propagation respectively at one extinction length from a 
300 nm laser source. As the detector approached the axis 
of propagation more photons were received. The photons 
along the phi prime equal to zero axis represent single 
scattered radiation. Many of the photons received were, 
however, multiply scattered. | 

Figures 4 and 5 are similar to Figures 6 and 7 only 
in the latter figures the receiver is located five extinc- 
tion lengths from the source. The percentages of multiply 
scattered photons increases greatly as the receiver is 
moved either off from the axis of propagation or away from 
the source. At five extinction lengths, virtually none 
of the radiation received is single scattered. This is a 
strong indication of the hazards of single scattering models 
in the middle uv. 

Figures 8, 9, and 10 represent the contours of equal 
photon flux for 250, 280, and 300 nm wavelengths. The 


label for each contour is 
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270 
RECEIVER LOCATED ONE EXTINCTION LENGTH 
FROM THE SOURCE. 0 = О TO 9 DEGREES 
A * 300 NM. FIELD OF VIEW = I80 DEGREES 


FIGURE 4. Photon Scattering Diagram (8' vs $') as 
Viewed by a Receiver Located 0 Degrees 
from the Axis of Propagation. 
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FIGURE 5. Photon Scattering Diagram (0' vs $') as 
Viewed by a Receiver Located 0 Degrees 
from the Axis of Propagation. 
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270 
RECEIVER LOCATED FIVE EXTINCTION LENGTHS 
FROM THE SOURCE. 0 =O TO 9 DEGREES 
A - 500 NM. FIELD OF VIEW = I80 DEGREES 


FIGURE 6. Photon Scattering Diagram (8' vs ¢') as 
Viewed by a Receiver Located 6 Degrees from 
the Axis of Propagation 
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270 
RECEIVER LOCATED FIVE EXTINCTION LENGTHS 
FROM THE SOURCE. @= 36 TO 45 DEGREES 
\ = 3500 NM. FIELD OF VIEW = 180 DEGREES 


FIGURE 7. Photon Scattering Diagram (0' vs $') as Viewed 
by a Receiver Located 0 Degrees from the Axis 
of Propagation 
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Ν 


Relative Flux = - озу]; 2TR (cos04 _ < 


(4) 


where 2 TR“ (cos9, - COS Ө.) is the area of that portion of 

a hemisphere located R extinction lengths fae the source 
subtended by 91 and 95; N is the number of photons pene- 
trating that area, and N. ls the number of photons leaving 

the source. Each successive contour of photon flux represents 
a decrease of one order of magnitude in the relative photon 
flux. Thus, for example, the contour labeled 2 gives the 
location at which the photon flux is ρα. photons per unit 


area, where unit area is one extinction length squared. 


E. FIELD OF VIEW CONSIDERATIONS 

Figures 11 and 12 represent the relative number of 
photons detected at a receiver located zero to nine and 36 
to 45 degrees off from the axis of propagation at one and 
five extinction lengths respectively from a 250 nm laser 
source. A sixty degree (full cone angle) field of view 
detector would receive seventy-five to eighty percent of 
the incident radiation for the 94 (zero to nine degrees) 
and the θ2 (36 to 45 degrees) cases at one extinction length. 
At five extinction lengths these percentages tend to remain 
stable. Thus, given a minimum required photon flux at a 
detector, it is possible to select the proper field 


of view for that detector. Conversely, utilizing a fixed 
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field of view detector, it is possible to determine the 
maximum effective range for a minimum photon flux. (The 
results for 280 and 300 nm were substantially the same 


as the 250 nm simulation). 


F. PULSE SPREADING CONSIDERATIONS 

Figures 13 and 14 represent the pulse spreading as 
received by a detector located zero to nine and 36 to 45 
degrees off from the axis of propagation at one and five 
extinction lengths from a 250 nm laser source, respectively. 
The pulse spreading function is a function of the position 
of the observer and time dispersion caused by multiple 
scattering. The relative time dispersion, T(R,D) is 


defined as follows: 


where D is the total distance traveled by each photon and 
R is the distance from the source to the shell penetrated 
by each photon. 

The results of pulse spreading at other distances 
and angles was examined. The pulse spreading 
function is nearly constant for values of R from one to ten 
extinction lengths, thus one figure suffices as an approxi- 
mation to any distance in this angular range. The increase 
in pulse spreading with increasing angle is illustrated 


with one additional plot at 36 to 45 degrees off from the 
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axis of propagation. Reproduction of the results at each 
cone angle introduces a cumbersome presentation of this 
information, hence only one representative sample is 


included. 


G. SUMMARY 

The absorption and scattering coefficients as well as 
the phase functions from previous experiments were uti- 
lized and the behavior of multiple scattered radiation was 
simulated using the Monte Carlo method. This model was 
used to study the characteristics of the source beam at a 
receiver located an arbitrary angle off from the axis of 
propagation at a given distance from the source. With such 
information available, the geometric pattern of possible 
communication links can be determined. A photon density as 
a function of receiver field of view and position of the 
observer was determined as well as the effects of multiple 


Scattering on pulse spreading. 
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VI. SIGNAL-TO-NOISE CONSIDERATIONS 


A. INTRODUCTION 

The photon flux information presented in the previous 
section may be utilized to estimate the signal-to-noise 
ratio (S/N or SNR). This information, however, must be 
used in conjunction with the characteristics of the detec- 
tor. Additional considerations therefore include the 
characteristics of the optical filter, quantum efficiency 
of the detector, background level of radiation, etc. This 
section considers the characteristics of available devices 
and of known atmospheric flux levels in order to estimate 


the flux level necessary for reliable communication. 


B. FUNDAMENTAL RELATIONSHIPS 
The detector current generated by a light signal of 


power P. and of frequency v is 


S 


| GneP. 
S hv 


where G is the gain of the detector, n is the quantum 
efficiency of the detector, e is the electronic charge, and 
h is Planck's constant [18]. 

The dominant noise mechanism for detectors operating 
in the visible and ultraviolet region of the spectrum is 
shot noise. The mean square of the noise current due to 


amis noise mechanism is 
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| d dues ne 
«i > = 2eG [ii E hv (Ps + Pp) ] B (7) 


where G is the gain of the detection device (typically a 


photomultiplier), i, is the dark current, P, is the power 


d B 
of the background radiation present due to sources other 
than that from the desired signal, and B is the bandwidth 
of the receiver [18, 25, 26, 27]. 


The desired expression for S/N is thus the ratio of the 


square of Equations (6) andof (7). 





ΠΕΕ 2 
p E 
S/N = --- = = гоган 505 (8) 
<i, > 28614 -- ο. (Pç + Pp) ІВ 


In several specific cases the dark current i, is negligible, 


d 


whereupon Equation (8) reduces to Equation (9). 


νη»... (5) 





For pulsed signals having a time duration of T, the band- 
width B is effectively 1/2т. The relationship for S/N 
thus reduces to the number of photons converted to electrons 
in a pulse, degraded by the ratio of the signal power to 
total power received over the duration of the pulse. 

Each of the expressions for S/N may be modified slightly 
by consideration of a duty factor, introduced to account 


for the fact that a signal is not constant over the duration 


48 


of a pulse. Consideration of this introduces a multi- 


plicative factor equal to Equation (10). 


Pis = 2 
μποτ 4 
05 _ -т/2 15 ο) (10) 
Doe 1/2 TE (СЕ 
2 f ——— dt 
-1/2 EE 


where De/Dy is the ratio of the duty factor of the signal 


to that of the noise over the duration of the pulse, τ. 


C. BACKGROUND LEVELS OF RADIATION 

The expressions for S/N involve critically the back- 
ground levels of radiation due to other sources. In communi- 
cation experiments this is typically due to scattered sun- 
light. Reliable flux levels at wavelengths longer than 
300 nm are reported by Valley [28]. The flux levels at 
shorter wavelengths are a result of both direct and diffuse 
transmission. These levels are critically dependent on 
the thickness of the ozone layer, and have been calculated 
for various atmospheric conditions [12]. 

The photon flux calculated by Shettle and Green is 
presented in Table IV. The dramatic drop in background 
radiation levels DO 280 nm is evident. The signal 
that would be transmitted through a narrow band (10 nm) 
multilayer dielectric filter at 300 nm would include up to 


ο Dhotons/cm^ as background noise photons, whereas at 
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wavelengths below 280 nm this background radiation level 


may be considered to be quite small. 


TABLE IV 
Photon Flux Due to Direct and Diffuse 
Transmittance through the Earth's 
Atmosphere [12] 


A solar angle of ninety degrees is assumed, along with a 


0.32 cm ozone thickness. 


À (nm) Flux απ σα nm) Flux (noon fame nm sec) 
300 1.07 x 10 * 1.6 κ 1012 
295 5.15 x 10 ^ 7.8 х 1010 
-6 8 
290 2.14 X 10 3.2 X 10 
285 9.50 x 10 11 1.4 х 10° 
280 9.09 x 10 12 ο... 


D. DETECTOR CHARACTERISTICS 

The detector characteristics which relate to a given 
S/N expression include filter characteristics, quantum 
efficiency, detector area, and field of view. Using a 
combination of multilayer dielectric filters and absorption 
filters, it ls possible to construct a device having com- 
posite characteristics which will pass ten percent of the 
incident flux at 265 nm and reject the ambient solar flux 
in the visible to a degree sufficient to reduce the overall 


background count level to approximately ten counts per 
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second [29]. The transmission characteristics of this 
filter in the 200 to 400 nm region are indicated in Figure 
15. It is thus expected that somewhat greater than one 
percent but less than ten percent of the photons in the 
250 to 280 nm region will pass through the filter. Once 
the photons reach the photocathode of a typical photomulti- 
plier, the conversion efficiency is usually greater than 
ten percent but less than 25 percent. Thus, for purposes 
of estimating S/N, the combined effects of filter trans- 
mission characteristics plus quantum efficiency of the 
photocathode may be accounted for by choosing a value of 
n to be greater than 1057 but less than i. Wlth unfor- 
seen developments it is possible to expect a combination 
having a composite value of ος however such a device is 
not presently available. 

The area of the detector used in the simulation was 


chosen to be nominally one square inch, typical of currently 


available photomultiplier tubes. 


E. STATISTICAL CONSIDERATIONS 
The number of photons converted into an electrical 
signal follows Poisson statistics: 


P(r) = ёз n"/r! (11) 


where P(r) is the probability that r photons are observed 
in a sample pulse which has an average number of n photons 


per pulse. 


SL 


TRANSMISSION 





200 300 400 
WAVELENGTH (NM) 


FIGURE 15. Transmission Characteristics for a Typical 
Composite Filter 
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The probability of observing an event with a number of 


photons above the threshold value Ton 15 


= 


) nt /r! (12) 


P(r?r 


um 2 2 
ΞΥ 


r=r th 


Graphs of this function are readily available [30]. Given 
an average number, n, equal to ten and a threshold level of 
two (to distinguish from the random single photon events 
which occur at the detector with a frequency of occurrence 
equal to ten or twenty photons per second), a probability of 
detection of such pulses is .9995. If the average number 

n is five with the threshold number equal to two, the 
probability of detection is 0.96. Thus an average value of 
ten photons per pulse in the absence of background radiation 
is chosen as a desired level to avoid excessive error. This 


means that a desirable value of S/N is ten. 


F. SUMMARY OF SIGNAL-TO-NOISE CONSIDERATIONS 
For the case of minimal background radiation, Eq. (13) 
eaves the relationship for the number of signal photons 


required. 
S "dcs 
mem (x! (13) 


If background radiation dominates the power incident on 
ehe detector (as would be the case at wavelengths longer 
than 300 nm), then the number of signal photons required is 


given by Eq. (14). 
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(14) 


Table V gives a summary of the necessary photon flux 


external to the detector in order to obtain a S/N equal 


to ten. 


These flux levels may be used in conjunction with 


the figures of Section V to generate the geometric patterns 


over which signal transmission may be expected. 


TABLE V 


Photon Flux Required in order to Give 


a S/N Equal to Ten 


Assumed parameters: 


Note: 


À (nm) 


250 


265 


280 


300 


Filter characteristics for Figure 15 


Quantum efficiency of ten percent 


τ 


=e ome 


seconds (significant pulse 


Area of detector is one square inch 


Filter characteristics at 300 nm have 


for a narrow band (ten nm) multilayer 





filter. 

e 

—B (photons, 
n hv sec 
107° --- 
ΠΡΙ. --- 
TONS --- 
José Ў 
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hv 


Pot 
(photons) 


— 
10 


103 


10“ 


102 


spreading) 


been assumed 


dielectric 


Required flux 
(photons/em^) 
1.6 X 10° 


1.6 X 103 


роб 10° 


1.6 x 104 


ντ. DISCUSSION 


Although not currently available, three millijoule 
laser sources of different wavelengths (250, 280, and 300 
nanometers) capable of being pulsed at an appropriate rate 
(200 pulses per second for 2400 bit per second data rate 
using 12 bits per pulse [4]) were assumed. (The navy has 
used vocoders with 2400 bits per second capability since the 
1940's). Table V lists the assumptions leading to the 
required photon flux in order to achieve a S/N equal to 
ten. Ву Utilizing the contours of equal photon flux in 
Figures 8, 9, and 10, the ranges for which pulse position 
modulation communication is possible were predicted. 

Table VI lists the results of these calculations. 

Predicted performance for lasers operating at 250 and 
280 nm are comparable for sixty degree or greater off axis 
detection. For detection angles less than sixty degrees 
the 280 nm laser is superior. The performance of the 300 
nm laser does not appear to be degraded by the background 
radiation as much as might be expected. The performance 
for this laser from zero to twenty degrees off from the 
axis of propagation is similar to the 280 nm laser. Beyond 
thirty degrees the ranges decrease markedly, indicating 
very little backscattering. It cannot be overstressed 
that the results above are for a simulation of assumed 


laser sources with assumed atmospheric propagation 
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characteristics and assumed detector characteristics. 
Clearly, in order to make accurate predictions for real 
systems it is necessary to obtain accurate atmospheric 


characteristics as well as other system parameters. 
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VIII. CONCLUSION 


Lasers operating in the middle uv do not appear to be 
suitable for long range communications; however, they pre- 
sent definite possibilities for short range (several km) 
applications. Since radiation in this region is absorbed 
exponentially, rather than as the inverse of distance 
squared, the potential for covertness is highly accent- 
uated. Current lasers, however, having a power output on 
the order of 0.1 joules or greater with low pulse rates 
(15 to 20 pulses per second), are not suitable for voice 
communication due to the low repetition rates. The pulsed 
mode of operation appears to hold definite advantages over 
amplitude modulation of a continuous wave laser from simple 
signal-to-noise considerations. Good possibilities for 
successful communication are offered by pulse position 
modulation. Millijoule laser sources operating with 200 
or more pulses per second are desirable but not currently 
available. 

The contours of equal photon flux presented in this 
paper are useful for determining the communication patterns 
of possible laser sources. The pulse spreading and field 
of view figures are likewise useful for predicting the 
characteristics of proposed receivers. A multiple scatter- 
ing model is definitely required for proper prediction of 


photon flux over expected distances of communication. 
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APPENDIX A 
MONTE CARLO SIMULATION OF THE 
MULTIPLE SCATTERING PROBLEM 


A. FUNDAMENTAL RELATIONSHIPS 
1. Weighting Factors for Monte Carlo Calculations 

The following calculations served as a model for the 
derivation of the various weighting factors: 

PROBLEM: Given a random number generator that 
provides numbers within a specified interval along the 
X-axis (with equal probability for each interval dx) 
obtain random numbers with a probability p(y). I.e., 

Find y(x) such that p(y)dy = dx. 
SOLUTION: Let Random Variable Y be distributed in 


DEUS Y < Y, and its density function be given by 


p(y) M Sa SY.) 
p(y) = 


0 otherwise 


πρ απ Y 1IDu ion ΕΙπσ πφη, Ρο), is given by 


O ` 


y 
P(y) = P(Y <y) = ο ο Gay 
y 
Therefore, P(y) = I iy) ay 62 р(у) = 0 ғог у < у 
Yo 
Now, 
ply) dy = ах 
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Integrating, 
— 1 
y(x) Sy EUM) 


where 


he 
II 


у(х.) 
апа Ax == Xo 
This concept applied to specific cases yields the 
following fundamental relationships for weighting factors: 
a. Exponential 


p(y) = p 


ar 


y (x) ^-T1 ln (1 - Ax) 


b. HenyeyzGreensteın 


(1-92) 
p(0) = INL (ο «ο Rn) 
2(ltg -2g cos 9 


where 


«cos Q0» 


Q 
lI 


8 (x) 


II 
Q 
O 
(n 


1+92 _ı[| TP 
29 2a 2g itg = | 
g 
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c. Modified Henyey-Greenstein [Zachor 1977] 


" E l DEMOS cos^ 0-1) 
NI Ú“ ποιο τσ 372° 
(1+g -2g cos о) (1+g ) 
О < оис T 
where 
ΓΟ and f is a weighting factor. 
2 2 2 
ΠΝ ο 159 _ |09<-1) 1 a | 
8g [5c Eu oy pon 
where 
Πε 3 
πιο} = a NU (cos 0 - cos' 9) 
(1+9) 
d. Uniform 
= κ: Quo =< т 
2π NE e 
ф(х) = 2T Ax 
e. Rayleigh 
3 2 
p(8) = g( + cos’ 6) ESSET 
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Ө(х) = cos ` ΠΕΡῚ 


where 
3 2 
" b b 
A = 2 + E + 1 
b b 
Να Ἠ E 
and 
b = ϐ Ax - 4 
2. Position and Coordinates in Photon-Fixed Coordinate 
system 
a. Axis Rotation 
Axis rotation is accomplished by the following 
transformation: 
x' cos A0 cos Ad —cos A0 sin Дф -sin ЛӨ x 
y' - sin Ad cos Ab 0 y 
z ' sin A0 cos A$ -sin A0 sin AQ cos A9 Ζ 
b. Axis Translation 
Axis translation is accomplished by the following 
transformation: 
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3. Angles (0,06,0',7£52') for a Photon Passing Through 
a Shelll 


A laser is located at the origin of the E jd k. 


coordinate system. The initial direction of propagation of 
the collimated beam is along the -ko axis. 0 is the angle 
between the position vector R (pointing from the location 
where the photon penetrates a shell toward the laser source) 
and the vector -k.. φ 15 the angle between the projection 
of -R onto the d. c plane (denoted by R and the 


^ 


vector i . 


О 
A photon is located at the origin of the Nn т 
V, and the UL, Ds p coordinate systems. As the photon 


passes through the shell, 0' specifies the angular direction 


of photon travel (along U,) relative to R (the line-of-sight 


position vector). .$' is the angle between the projection 
of 0, (denoted by the vector W) onto the РОЈ, plane and 
NN. 

x 


The vectors R;, Ris and R. are position vectors 


in the Ur ο. U, coordinate system denoting the location 


^ 


of the points (1,0,0), (0,1,0), and (0,0,1) in the іс’ Jo 


1 сее Figure 16. 
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FIGURE 16. Diagram of the Photon Fixed Coordinate 
System 
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^ 


k. coordinate system. (These points are used for orientation 


purposes.) " д 
| T ) 
-1 O 
a. 0 = cos —— n 
IR| [kol 
IR (R. - R) 
8 = οκ | = (Om ο Ἢν 
к. і 
р ф = sos™? ( 2) 
R 
p 
[E Sw 
ф = cos = 
where 
- EN ο ο (0. 6 « 2m) 
το Ξ I K S Es 
= En 
С Ө = = 
Aa) 


ο o. scm) 


г) 
=, 


where 


4. Mapping SphericalCoordinates onto Polar Coordinates 
and Coordinate Transformation 


In order to map a point on a sphere onto a circle(s) 
such that there is a one-to-one mapping the following 
transformations were utilized: 


Case I. 0 < 0 < T⁄2 


r = R Vl - cos Ө 


Case [I πι or 


b 
II 
"Jj 
= 
+ 
Q 
O 
(n 
CD 


cD 
! 
Q 
О 
(n 
I 
= 
u 
AIR 
N 
i 
к— 
x) 


In order to plot (90,96) the following coordinate 


transformation was used: 


I.e. - σἰνθη (θ,ϕ), find (x,y). 
la - cos 6 0 < 0 < mi 
ү == 
R yl + cos Ө ПК о 


νά 
il 


r cos » 


r sin ф 


he 
N 
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5. Phase Functions 
A factor of i arises from the normalization 
process for the phase function. 
2T 


T 
(i.e. Ј πιο sintiderad = 19 
0 0 


a. Particulate 


B EG l 
P) = An | VU I-A EPI 


IFG -2g cos s 





ЕЗ СО 
PADRE EJUS 
2 TS ) 2 2 


4 


b. Rayleigh 
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APPENDIX B 


DEVELOPMENT OF THE FUNDAMENTAL RELATIONSHIP 
FOR RAYLEIGH SCATTERING 


p(9) = š (1 + σοβ΄θ) ος θη 


θ 
f p(90) d cos Ө 
0 


er 
ii 


9 3 2 
- Ј g(t * cos"0) d cos Ө 
0 
3 θ 
ath =) eens none Oj 
8 3 0 
3 ο. 9 T! 
x = cos ο μα 1 
Now solve for cos Ө. 
Re, + 3 cos 6 + (8x - 4) = O 
Let 
y = cos Ө 
b = 8x- 4 
a = 3 


68 


u i >a 


(See Reference 3l for the solution to a cubic equation.) 











у = A+B 
ог 
ро = - зв ииз "7ч 
Or 
where 
B - Vey аз 
and 
aN ια мш 
2 4 27 
Now 
S e 
4 27 4 2η 4 
ο΄ EO be 'ὗ 10 


Therefore, Only one real root and two conjugate imaginary 


roots exist. 


Since y = cos 9, only one real root is considered. 


Coe ο + B) 


o = 
where 
3 2 
_ b b 
3 2 
B b b 
PIS NI re 
and 
b = 8x - 4 
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APPENDIX € 


SAMPLE CALCULATION FOR RANGE ΡΕΕΡΤΟΤΤΟΝ 


PROBLEM: Given a single channel 280 nm laser utilizing 
pulse position with power output per pulse equal to one 
millijoule, determine the communication pattern for a S/N 
equal to ten. Assume the detector characteristics presented 


in Section VI. 


SOLUTION: From Table V, the minimum photon flux for a S/N 


of ten is 1.6 x To photons per em”. 


Now, 
] watt = 5.0345 X 1018 λ photons 
sec 
where 
A= 2280 1 (A is meastired in u) 


1 joule = 1 watt sec 


1.40966 x ο photons 


1.40966 x ο photons 


l millijoule 
From Figure 9, at contour 2 the relative fluxis Ге of the initial 


photons per unit area where the unit area is one extinction 


length squared. From Table TII one extinction length is 
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(1.0405 km 1j 1 or .9611 kilometers. The relative flux is 


1.0826 x ο. cm ^. The irradiance is given by the product 


of 1.40966 x 107 апа 08 E US 


or 1526 photons per 
Square cm. By measuring the distance from the source to 
the detector located on contour 2 (Figure 9) and recalling 
that one extinction length for this case is .9611 km, the 


ranges for an irradiance of 1526 photons/cm“ may be 


determined. Table VI lists the results of this calculation. 
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APPENDIX D 


CHECKS ON POSSIBLE ERRORS 


Each computational procedure involved in writing this 
program was verified and an effort was made to ensure the 
entire program functioned properly. The procedures uti- 
lized to document the correct behavior of this routine are 
described in this appendix. 

The various weighting functions in the program were 
verified by comparing the actual distribution of photons 
generated by the computer to the analytical solution for 
the distribution of random numbers obtained using a program- 
mable calculator. (The statements required for this genera- 
tion remain in the program, but the output of the results 
is suppressed.) 

Verification of axis transformation computations was 
simple and straight-forward. The axes were rotated about 
the Z axis, the Y axis, and finally about both axes to 
ensure proper behavior. (No documentation steps of this 
nature have been left in the program.) The output of 6, 
$,0', and $' for each photon at each shell was ordered by a 
computer sort routine in order to facilitate the preparation 
Of a scatter diagram indicating the pattern formed by the 
simulated photons. Inspection of such a pattern did in 


fact reveal the presence of a fault which was corrected. 
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The results of Duntley [32] in describing multiply 
scattered radiation in the propagation of laser light under 
water are most useful. Flux as a function of cone angle is 
presented for distances of up to 19 attenuation lengths. 
The scattering phase function for light underwater is a 
sharply peaked function which may be approximated with a 
Henyey-Greenstein function with a g value equal to 0.96 
[33]. Calculations of photon flux as a function of cone 
angle were carried out with the computer program developed. 
Total attenuation of the beam down to a level of ΙΩΝ. could 
be calculated with some measure of accuracy (only 10,000 
simulated photons were utilized for these calculations.) 
The results of these calculations agreed with the behavior 
observed in these underwater scattering measurements, simu- 
lating correctly the flux levels out to 16 attenuation 
lengths (the limit calculated). 

Further verification of the accuracy of the program 
was obtained by comparison with calculations by Zachor and 
Green [24]. In this article, photon flux as a function of 
viewing angle was calculated for several distances from the 
source. The same conditions were assumed with the model 
used here, resulting in satisfactory agreement with the 
information published by Zachor for distances of 1.43, 2.73, 


and 4.95 extinction lengths for 300 nm radiation. 
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APPENDIX E 


PROGRAM DESCRIPTION AND DOCUMENTATION 


A. COMPONENTS OF THE PROGRAM 
J —Wesghting s Funcezon Subprogmsams 

RANEXP is a function subprogram that generates 
random numbers weighted exponentially. It is used for 
distance calculations. 

Function subprogram RANTH generates random values 
for theta and theta prime. Since the interaction may be 
either Rayleigh scattering or particulate scattering, 
either calculation may be performed. A uniform random 
number (generated by RANDU) within the interval [0,1] is 


obtained and compared to R_ (the ratio of particulate to 


Р 
total scattering). If the random number is less than R5 
the particulate scattering calculation is made. (Otherwise 
the Rayleigh calculation is performed.) Due to the com- 


plexity of the particulate phase function it is difficult 
to obtain an exact solution; therefore, an iterative method 
is employed. An exact solution is available when f is set 
to zero and a more efficient calculation is utilized. 

RANPH generates a random value of phi or phi prime 
uniformly weighted on the interval [0,27]. The steps in 


this function subprogram are trivial. 
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2. LITE Subroutine 

The main subroutine of the Monte Carlo simulation, 
called LITE, simulates a photon which randomly interacts 
with atmospheric constituents. (In order to make the most 
efficient use of computer time, non-interaction was not 
considered.) Initially the photon travels along the -Ko 
axis (see Figure 16, Appendix A), a distance AR (determined 
by function subprogram RANEXP). A uniform random variable 


in the interval [0,1] is compared to R. (the ratio of the 


S 


total scattering coefficient, K to the total extinction 


S” 
coefficient, K). If the random number is smaller than 


R the photon is scattered. (Otherwise it is absorbed and 


с” 
the photon is lost.) If it is scattered, the distance of the 
photon from the laser is updated and a new value of R is 
calculated. Since the photon has been scattered, new angles 
and distances must be computed. This is accomplished by 
generating values for Ad, Ad, and AR (using RANTH, RANPH, 

and RANEXP, respectively). The photon coordinate system is 
transformed (translated and rotated) to a new coordinate 
system. The distance from the source is calculated and 
compared with the radius of a shell. (A series of concen- 
tric shells with radii measured in extinction lengths, may 

be formed in order to determine the photon flux or angular 
distribution of photons as a function of distance from the 
source.) Various calculations are then made to determine 
whether a photon has penetrated through a shell, passed in 
and out of a shell, or penetrated an outer shell. These 


calculations are necessary to ensure proper accounting of 
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the photon flux. When a photon passes through a shell, the 
position of the source with respect to the penetration 
point is obtained. Theta, phi, and the radius of the 
shell (measured in extinction lengths) are stored in a 
three dimensional array for future use. Theta prime and 
phi prime are also determined as a function of distance 
(using RANTH and RANPH). These too, are stored in a three 
dimensional array and specify the direction of incident 
photon flux as observed by a receiver located a given dis- 
tance from the source at a given angle off from the axis of 
propagation. This process is repeated until all of the 
photons have been used to generate information. (It should 
be noted that many extra steps incorporated to debug the 
program, remain in to facilitate future implementation by 
another user.) 
3. DRLITE Routine 

The routine DRLITE is used to drive the subroutine 
LITE. Information is read into the computer and the results 
are written out using this routine. The spatial distribu- 
tion of photons as a function of position and detector field 
of view iS calculated as well as the total number of photons 
in each shell. The output of other available information is 
optional and is controlled by IPRT and METHOD statements. 
This information includes pulse spreading (or photon time 
of arrival) as a function of position and field of view, 


the relative photon flux at each annular ring ina shell, 


τὴ 


and a least square curve fitting of a modified Henyey-Greenstein 
function to the phase function information generated by the 
computer. 

4. Miscellaneous Subroutines and Function Subprograms 

LSTHG, LEAST, GAUSS3, and RANDU are miscellaneous 

subroutines used in the program.  LSTHG is used to least 
squares fit a modified Henyey-Greenstein function to the 
information generated. LEAST and GAUSS3 are used in con- 
junction with LSTHG.  EQN is a function subprogram supporting 
LSTHG.  RANDU is a uniform random number generator that 
provides numbers in the interval [0,1]. It is located in 
the computer library. The variable IX is the seed and may 
be any odd integer with nine digits or less. It may not 
be zero. In order to use RANDU the following FORTRAN 


Statements are used: 


IX = 948752759 
CALL RANDU (IX,IY,Y) 


IX = IY. 


The seed IX is used in the calculation and a new seed IY is 
generated for future use. A uniform random number, Y, is 
in the interval [0,1] and is also available for future 
calculations. The last statement (IX = IY) resets the seed 


IX to a new value. 
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